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1. Introduction 



Since tlie classic work of Candelas, Horowitz, Strominger and Witten in the mid 1980s 
1^, Calabi-Yau compactification of lOd supergravity has been the dominant paradigm for 
the construction of 4d string models. In the heterotic string theory, such compactifications 
yield M = 1 supersymmetric theories below the KK scale, and can quite plausibly acco- 
modate the standard model. In the type II context, Calabi-Yau compactification yields 
M = 2 supersymmetric models, which have been of considerable interest in studies of string 
duality and mirror symmetry. 

More recently, it has become clear that these models are only the simplest represen- 
tatives of a much larger class of models. For concreteness, we will restrict our discussion 
to the type IIB theory, but similar remarks would apply in the other limits of M-theory. 
In the IIB context, the theory contains p-form gauge fields which couple to the various 
supersymmetric D-branes and NS-branes. One can consider magnetic flirxes of these gauge 
fields in the internal dimensions; such fiuxes are consistent with 4d Lorentz invariance. In 
fact, in type II orientifold models which break the supersymmetry from A/" = 2 to A/" = 1, 
generic methods of satisfying the tadpole conditions will involve the addition of back- 
ground magnetic fiux. This is interesting because the magnetic fiuxes will cost an energy 
which depends on the detailed shape and size of the internal dimensions; that is, they will 
generate a potential for some of the Calabi-Yau moduli. 

In the case of type IIB on a Calabi-Yau orientifold, the allowed fiuxes are those of 
the RR and NS three- form field strengths F3 and H^. The properties of the theory in 
the presence of such fiuxes have been investigated in a vast literature; see for instance 
|^-|^. The resulting M = 1 supergravity theories (derived in [^,^) are of the no-scale type, 
with a Kahler potential K inherited (at large volume and weak coupling) from the Calabi- 
Yau compactification and a superpotential which is computable in terms of the choice 
of quantized background magnetic fiuxes 0. The resulting scalar potential generically 
fixes the complex structure moduli and the dilaton. It is possible to make semi-realistic 
chiral models of particle physics in this framework , and to describe some aspects of the 
resulting supergravity theories in the general framework of gauged extended supergravity 



||10|| . The theory of soft susy breaking terms in these models has been developed in . 

For a given choice of the integral background fiuxes, the superpotential for moduli 
can simply be written as a combination of periods of the holomorphic three-form O on 
the Calabi-Yau space. The periods are computable by standard techniques in classical 
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geometry, highly developed in the physics literature after the seminal work of Candelas et 
al on mirror symmetry [|12[. The resulting superpotential is exact at least to all orders in 
the a.' expansion by standard holomorphy arguments. It can be argued that "generic" non- 
perturbative effects (known to occur in various models in various circumstances, and likely 
to occur in a generic setting due to the absence of symmetries that would prevent them) 



should suffice to fix the remaining (Kahler) moduli in some cases [jlSiP^ . As argued in ||Ti 
if the SUSY-breaking (e^|VFp evaluated in vacuum) by the fluxes is suitably small, and 
if is stabilized by the fluxes at weak coupling, one can have small expansion parameters 
which allow one to stabilize the remaining moduli in a regime of control. Of course, this 
will only happen in at most a small fraction of the vacua. As emphasized in , following 
the basic ideas of Bousso and Polchinski [jl5[(see also [jl^), the large number of possible 
integral flux choices makes it plausible that suitable flux choices do exist, which allow one 
to flx at weak coupling and keep the SUSY breaking hierarchically smaller than the 
string scale. 

This set of issues, which goes into establishing the existence of a "landscape" of 
string/M-theory vacua (like the "discretuum" of [|1^) and characterizing its properties, 
can and should be studied systematically and quantitatively. In addition to the works 
already described, a persuasive discussion of this can be found in the works of Douglas and 
collaborators [p^p!7| -|T9| (see also pO| , |2T[] for more general comments about the landscape) . 



In particular, in [^, under well motivated assumptions some theorems are proven regard- 
ing the numbers and distribution of flux vacua on the Calabi-Yau moduli space (really, the 
moduli space of complex structures and dilaton VEVs) M.. The assumption which goes 
into proving these theorems is, roughly, that the superpotential is a random section (with 
respect to a covariance kernel determined by the Kahler potential) of the appropriate line 
bundle over M.. This models the flux superpotentials as random combinations of periods 
with respect to some measure, instead of as the precise integral combinations of periods 
determined by various choices of ifs and -F3 consistent with the tadpole conditions. Ba- 
sically, this assumption follows from the statement that if one allows large enough fluxes, 
flux quantization becomes irrelevant for the average properties of the vacua. The natural 
size of the fluxes required by tadpole conditions in IIB string theory probably falls into 
this range. 

Because the emerging discretuum of vacua is too complicated to allow us to usefully 
examine each vacuum individually, this kind of program to find simple assumptions which 
accurately characterize the space of vacua seems to us rather important. In fact, following 
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|T^,|T9[, one could simply disregard the Kahler moduli altogether, ignore the no-scale struc- 
ture (which is certainly broken by quantum corrections) , and characterize the properties of 
the supergravities that result with the given K and as a function of complex structure 
moduli. The resulting ensemble of supergravity theories is rich enough that it is very plau- 
sible that some aspects of the full story (about the real space of string vacua) , are captured 
by this simplified ensemble. However, we will continue to use the language and equations 
appropriate to thinking of these as vacua which arise in the no-scale approximation to the 
full IIB theory, though it is easy to translate back and forth. 

In this paper, we test the prediction of regarding the distribution of vacua on 
M., by studying the vacua which result on a subspace of the moduli space of a certain 
Calabi-Yau threefold, in the presence of a small subset of the possible fiuxes. While we are 
looking at a very small fraction of the possible fiuxes and hence of the total number of fiux 
vacua on this space, the numbers even on our subspace are easily large enough to allow a 
meaningful test of the conjecture. In addition, in these vacua there is a tadpole condition 
that requires / A F3 = A^Aux < L, where L is some number fixed by global details of 
the model. The growth of the number of vacua with L, fixing the number of cycles which 



support fiux, is subject to various conjectures along the lines of [|T5|. Here, in addition to 
testing the distribution of vacua on moduli space, we also test the growth of the number 
of vacua as a function of L, and find a simple power law. A proper understanding of the 
theory which determines this power is rather complicated, as we describe in the following; 



some relevant results are discussed in the revised version of fl9 



The organization of this paper is as follows. In §2, we describe the predictions we wish 
to test. In §3, we discuss the model in which we will test them (which is "model A" of 0). 
This model has (on the locus of interest to us) a one-dimensional complex structure moduli 
space, parametrized by a complex variable tj). In sections 4 and 5, we describe how we find 
solutions of the fiux equations in two regions in the moduli space, using expansions of the 
periods around '0 = (the mirror Landau- Ginzburg point) and -0 = 1 (the conifold point). 
In each of these cases, we find sufficiently large numbers of vacua to perform meaningful 
tests of the predictions. In §6 we present our conclusions. 

We should stress that although the discussion here is specific to type IIB vacua of the 
class described in Q (for which several de Sitter scenarios exist [jl^,^), a similar landscape 
also emerges from careful consideration of the various other limits of M-theory (see e.g. 



for discussions in the heterotic setting, and ||2^ for an early construction of de Sitter vacua 
in noncritical string theories). For instance in the heterotic string, one needs to consider 
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the large class of possible gauge bundles in each Eg, in addition to the possibility of NS 
three-form flux and Non-Kahlerity . This is technically more challenging than studying 
the class of IIB flux vacua we considered here, but flnding similar tractable toy ensembles 
in that setting should be possible. Another setting where similar equations arise is in the 
problem of flnding BPS black holes in A/" = 2 supergravity; the interesting connections 
between that problem and the problem of flnding flux vacua are described in . 



While this work was being completed, the paper |19] which has some overlap with our 
results appeared. 



2. The conjectures 

We will be interested in testing two conjectures about the large class of IIB flux vacua 
described in [^, in the no-scale approximation. The conjectures concern the distribution 
of flux vacua on the moduli space of complex structures, and the scaling of the number 
of flux vacua with the D3-charge in the fluxes. They can be proven to hold [jl^ under 
suitable assumptions about the irrelevance of flux quantization, so we are really testing 
whether those assumptions hold for simple representatives of the class of vacua described 
ini. 

2.1. Where do the vacua lie? 

In a recent paper of Ashok and Douglas [jl^ , it was conjectured that the index density 
of the no-scale vacua described in is given by an integral 

J~ / det(-i?-w) . (2.1) 
Jm 

X is normalized appropriately to the total number (or index) of no-scale vacua. Here uj 
is the Kahler form on the moduli space and R is the curvature two-form expressed as a 
Hermitian matrix 

Ri = Rgi^dz' A dz^ (2.2) 

with the coordinates on Ai. The moduli space Ai consists of the fundamental domain of 
any relevant duality groups, acting on the complex structure and axio-dilaton Teichmuller 
spaces. In our example, it will be a product of the standard fundamental domain of 
SL{2, Z) in the upper half-plane (for the axio-dilaton r) with a (complex) one-dimensional 
slice of a Calabi-Yau complex structure moduli space. As described in [ir7| , p!9[] (see especially 



§3.1.2 of [|T^), one can integrate over r-space and find a distribution of vacua on just the 
complex structure moduli space; for one-modulus models one finds 

X~ / ci , (2.3) 

JM 

where now M. is just the complex structure moduli space, and ci is the first Chern class. 

We will assume that there are no miraculous cancellations, so that the total number 
of vacua scales like X. In this case, one predicts that the number of vacua in a given region 
of complex structure moduli space should scale as in (p.3|). We will test this conjecture 
in sections 4 and 5 by computing the numbers of vacua lying in various regions of the 
moduli space (by using monte carlo simulations to choose fluxes and then solving the flux 
equations), and checking that they are distributed like the appropriate integrals of the first 
Chern class. 

2. 2. Growth of number of vacua with L 

The three-form fluxes in type IIB supergravity contribute to the D3-brane charge 
tadpole on the compact space, via an equation of the form 

where tadpole cancellation requires 

iVD3 + iVfl,, = L . (2.5) 

Here L is some number fixed by the orientifold charge in the given model and A^ds is the 
number of space-time filling D3 branes, transverse to M. 

While it may appear from (|2.5D that this is not much of a restriction, since the 
quadratic form J A is not positive definite, this is not quite correct. On solu- 
tions of the equations of motion, one finds (as on p. 4 of the first reference in §) that 
-^3 ~ — *6 so Afflux becomes a positive semidefinite form. This means that in a space 
with K three-cycles, one might heuristically think of ( p.5|) as allowing one to search for flux 
vacua at each integral lattice point within a sphere of radius \/L in flux space. Although H 
and F together involve 2K flux integers, using the relation between H and F on solutions 
indicates that for a fixed point in moduli space, this will he a K — 1 sphere (instead of 
a 2K — 1 sphere). So for e.g. a (T^)^ with square two-tori, one can easily rewrite ( |2.5| ) 
(setting No3 = 0) as the equation for a sphere in K dimensions. Then, the number of 



5 



solutions for the flux equations Nyacua, like the number of lattice points inside this sphere, 
presumably scales like 

N.acua ~ L""/^ . (2.6) 

This would be the behavior of the number of solutions, in the simple toroidal example 
above, at a point in Ai where solutions to the flux equations exist. 

This is the kind of logic suggested in e.g. applied in this context by many 

workers [^, and rigorously argued to hold in this context in [|T^. However, the results of 



contain a surprise - the predicted scaling of the total number of vacua is actually 

N vacua ~ (2.7) 

instead of the L^/^ in ( |2.6| ). To understand this,l one must carefully consider a smoothed 
average of the number of solutions per unit volume in M. - the number of solutions jumps 
disco ntinuously as one moves in , so one must average to get a reasonable formula. This 
averaging produces an enhancement of the result - for large enough L, a given small region 
of moduli space will give rise to lattice points within flux space, instead of the naive 
fpj^jg jj^ particular means that if one looks at flnite L in some region of the moduli 
space, one may flnd an intermediate scaling that is between the naive prediction (|2.6|) and 



the behavior (p.7| ) expected from JT^. The detailed crossover behavior determining 



how the scaling should go for a given region in M. at various L and K is a. complicated 
subject which is still under investigation (see the revised version of the paper by Denef 



and Douglas [|T^ for further discussion). 

Because of the complications described above, we simply take our results as an exper- 
imental clue that should help in the development of the full theory of how Nyacua scales 
with L and though we expect to (and do) flnd a power intermediate between (|2.6| ) and 



^ One might ask whether one could increase the number of solutions by including anti-D3 
branes (i.e. taking Nds < in ( |2.5| )). While one can include anti-branes in these models, the 
number is bounded by the decay mechanism discussed in p7[] , and we will simply assume Nd?, > 
for now. 

^ We are grateful to M. Douglas for several discussions about this subject. 
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3. The model of interest 



We will study an orientifold of the Calabi-Yau threefold M which arises as a hyper- 
surface in the weighted projective space ^^-^1^,1,1,1,4 



4:Xq + Xi + X2 + + — 8tpXoXiX2X3X4^ = 0. (3.1) 

This threefold has /i^'^(M) = 1 and /i^'^(M) = 149. As explained in [0, an orientifold 
of M arises as a limit of the F-theory compactification on the Calabi-Yau fourfold X 
given by a hypersurface in W^-Pf,i, 1,1, 8,12- The orientifold action takes xq — > —xq and 
'i/' —>■ —-0 (in conjunction with the worldsheet parity reversal), while leaving all other 
variables unchanged. This limit of the F-theory model on X is obtained by following the 
general prescription given by Sen [^^. One of the benefits of knowing the F-theory lift 



is that it allows us to easily compute the tadpole condition for compactification on the 



orientifold of M. The relevant condition is (see and the third paper listed in [0) 

iVD3 + Nan. = ^ = 972 (3.2) 

where A^ds is the number of space-filling D3 branes one inserts transverse to M, and Afiux 
is the amount of D3 charge carried by the and F3 fluxes as given in (|2.4|) . 

In general, one could add many deformations to the deflning equation ( p.l|) which 
are consistent with the orientifold action. However, as explained in following |]3T[] , 



the equation ( p.l|) is the most general equation consistent with a global symmetry group 
Q = Zg X Z2 acting on the projective coordinates. If one turns on only fluxes which are 
consistent with preserving this symmetry, then the resulting superpotential is guaranteed 
to depend on the additional moduli only at higher orders (in Q invariant combinations). 
This means that if we set these moduli to zero (as done in (|3.1[) ), compute the periods 
relevant for our subset of fluxes, and solve the resulting equations for i/j and the axio- 
dilaton r, we are guaranteed that the solutions we flnd will be solutions of the full theory. 
The additional moduli can be consistently set to zero because of the symmetry, and will 
generically be constrained by a higher order potential in the resulting vacua. 

In fact, the subset of periods which are relevant to the Q invariant fluxes coincide (up 
to a factor of with those of the mirror manifold to M, W, as explained in 0. These 
correspond to periods on four different 3-cycles in M, and hence we proceed to develop 
the relevant formalism for computing flirx superpotentials involving only and F3 flirxes 



7 



through these four 3-cycles. Happily, the periods of interest to us in this example were 
computed in one of the early studies of mirror symmetry, by Klemm and Theisen 

In comparing to the conjecture of §2.2, we see that our problem has K = 4 and 
L = 972. So in the physical problem of interest, we could check the number of vacua with 
flux charge A^Aux < for any L* < 972. (Of course formally, we are free to consider vacua 
using the periods from this model with arbitrary L). 

3.1. Homology and Cohomology bases 

We will work with a symplectic homology basis for the subspaces of Hs of interest to 
us. The basis of three-cycles Aa and (a = 1, 2) and the basis for integral cohomology 
CKa and (3°' satisfy 

/ at = 5^, [ /3« = -5,^ / aaAf3' = 5l (3.3) 

JAa JBb JM 

The holomorphic three form can be represented in terms of periods in this basis as follows: 
In addition 

J M 

Here, we have introduced the prepotential Q{z^,z^), the period vector H (whose entries 
are the periods (|3.4|)), and the matrix 

whose entries are two by two matrices. 

3.2. Fluxes, Superpotential and Kdhler Potential 

The NS and RR fluxes admit the following quantization condition 

F(3) = {2'Kfa'{fa[Ba] + /a+2[A„]), if(3) = {2Tifa\K[Ba] + /l„+2[A„]) (3.7) 

with integer fi and hi. Here we used the notation [Aa\ = aa and [Ba\ = (3a. Using this 
notation, we flnd the following expression for A^Aux 
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O = Z'^aa - Gals'". 



(3.4) 



^«^_;,a^ = _nt.E-H. 

dz°' dz"" 



(3.5) 



The superpotential is given by 

W = [ (F(3) - rif(3)) A O = (27r) V(/ -U-rh-U). (3.9) 
The Kahler potential for the dilaton-axion and complex moduli is given by 

K = Kr+K^ = -ln{-i{T-T))~ln{i [ OaO) = - ln(-z(r-r)) -ln(-zn^ -E-H). (3.10) 

3.3. Conditions for solutions 

The tadpole condition for the model of interest is given by ( |3.2| ). This provides a 
restriction on the fluxes we can use, i.e. on our allowed class of superpotentials - we must 
keep Afflux < 972. 

We will be searching for nonsupersymmetric flux vacua which obey the equations 

D^W = and D^W = 0. (3.11) 

They are nonsupersymmetric because in the full IIB string theory, they will generically 
have DpW ~ 7^ where p is a Kahler modulus. Note however that in the toy ensemble 
of [lT7|,p!9[, such vacua are called supersymmetric, since the Kahler moduli are not part of 



the toy ensemble. 
4. Small i/; region 

4.1. Periods, Kahler potential, Metric and Curvature form 
The periods for the small '0, < 1 region are given by 

* ' Sf-.rWlm-DWl-i) sin i, 
where a = exp (^). The corresponding period vector is given by 

We will want the periods in a symplectic basis. In such a basis, the period vector Yl'^^i/j) = 
{Qi{4'),Q2{i^), z^ii^), z'^ii^)) is provided by the linear transformation 

n = m-w, (4.3) 
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where the matrix m is given by 

/ 2 2 2 2 \ 

0-10 
-1 3 2- 
V 1-10/ 

The period vector in the symplectic basis can be written in the form 

oo 

n(V;) = 5^C2nP2nV''"- (4-5) 



(4.4) 



1=0 



Here the vectors p2n are given by 

P2„= m-(a2(2n+l)^^(2n+l)^;L,«^(2n+l))T ^^q^ 



(4.7) 



and the first few constants C2n of interest are given by 

^° - ^'^^^ 2r4(7/8) sin(f) ' - -^'^^^ r4(5/8) sin(f ) ' 

C4 - (2.:)3 exp(3p) _ 3 exp(If) 

This yields the following Kahler potential for the complex structure moduli space: 

K^ = - ln[Ao + AalVI^ + + AsV'^ + + 0(|V;|^2)] ^4^8) 

where the real constants are given by 

Ao = -i|coppJ ■ S - Po, A2 = -i|c2pp5 ■ S ■p2, 
A4 = -^|C4| P4 • S • P4, A5 = -tcoCsPo • S • Po- 

Note that pj • E • po = 2i{V2 + 2) and • E • p2 = -pl • S • ^4 = 2i{V2 - 2). 
The resulting Kahler metric takes the form 

9^^ = 1^1^ (1 + (4^ - 2^) IV^r + 0(V'«)) . (4.10) 

Using the expression for the curvature form 

R^^ = Ki^^--9^^H9^^) (4-11) 

one finds the following expression of interest 

R^^=-M^\'(^¥-^¥)+Om% (4.12) 



V^- -1^1 V A^ A^ 

Note also that one finds 



d^K^ = -2^\ij\'',l; + 0{iP^), (4.13) 



which is of use in evaluating Kahler covariantized derivatives with respect to 
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4-2. Distribution of flux vacua 

Given these facts about the Calabi-Yau moduH space geometry, it is now possible to 
work out the detailed prediction from §2.1. One sees that 



2T:ici{ip) — const / {ipl' 

M Jo 



const - 



(4.14) 



for (an eight-fold cover of) a small piece of the fundamental region of complex structure 
moduli space M.: li/j] < r. 
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Fig. 1: A plot of a numerical evaluation of 2TTici{ip). The monte carlo sim- 
ulation data for each point is: number of random fluxes = 10*; random flux 
interval f,he (-40, 40). The data is fit by the curve 700000 ^. 



The equations for flux vacua, keeping terms up to order ij;^, become 

DrW = {f-Th)-{coPo + C2P2i^^) = 0; 



,2 ''^2 —2s 



D^W =0 =^ U-rh)- (C2P2 + 2C4P4V' - T^coPoV' ) = 0. 



This provides two equations 



—2 



h-pl + f^h-pl^ h.p2+^h-p^^^-^h-p^^ 



(4.15) 



(4.16) 



11 



The last one gives a linear equation for ip'^ such that 



a + bii;^ + baV'^ + 0(V'^) = 0, 



(4.17) 



where we defined the following quantities depending on fiiixes 



a = f ■ Poh ■ P2 - h ■ pI f ■ P2; 

bi = —{h ■ p4 f ■ pI - f ■ p^h ■ pI); 

69 = ^^{f ■ Po h ■ pI - h ■ po f ■ pI) + ^{f ■ plh ■ p2 - h ■ pI f ■ P2). 



(4.18) 



A0C2 



Co 



Introducing the notation i/j'^ = x + iy one gets the following system of linear equations 



Re(a) + Re(6i + 62)2; - Im(6i - 62)2/ = 0; 
Im(a) + Im(6i + b2)x + Re{bi - 62)?/ = 0. 



(4.19) 
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Fig. 2: The numerical results for the number of vacua with A^aux < L for L G 
(1,972). The monte carlo simulation data is: number of random fluxes N = 10®; 
random flux interval f,h E (—100,100); complex structure tp space region I^P < 

r 3 

0.5. The data is fit by the curve 



183000 • 
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We used monte carlo simulations to check the validity of the formula ( [4.14 ) and the 



scaling of the number of vacua with L. We chose correspondingly 10^ and 10^ random 
fluxes, taking integer values /, h in the ranges indicated in the figure captions. The number 
of solutions with L < 972 resulting in small is determined numerically from eqs. ( [4.16| )- 



( [4.19| ) and is plotted in Fig.l. Fig. 2 shows the dependence of the total number of vacua 
on L. 

One sees that the curve in Fig. 1 is in very good agreement with the prediction ( |2.1| ). 
It is interesting that in Fig. 2, one sees a rather clean growth of the number of vacua. 
The naive formula (|2.6|) would predict L^, while the more refined formula proved by Ashok 



and Douglas under certain assumptions would predict L^. In fact, as described in §2.2 and 



the revised version of [|T9[ , the correct scaling for a small region in -i/^-space is a complicated 
subject. It may be that exceedingly large L would be required to attain the scaling ( |2.7|) , 
in this region of moduli space. 

5. Conifold region 

5.1. Periods, Kdhler potential, Metric and Curvature form 

The periods in a symplectic basis in the vicinity of the conifold point = 1 can be 
given to the first order by the following expressions (here x = 1 — and |a;| ^ 1) 

Gi{x) = (27rz)^[ao + aix + 0(x^)], 

G2{x) = ^ Hx) + {2mf[bo + bix + 0{x% 
z^{x) = (27rz)^[co + cix + 0(x^)], 
z^{x) = {2nif[do + dix + 0{x'^)]. 

where the constants can be approximated by the following numbers 

ao = -1.774i, cq = 4.952 - 5.321z, 

ai = 1.227z, Cl = -4.488 + 3.682z, 

bo = -1.047, do = 0, 

bi = 0.451 + 0.900Z, di = 1.800i. 

The Kahler potential for complex structure modulus is given by 
= — ln[no + Hix + JIiX + i-i2\xf In |xp + /^ala^l^ + /i4a^^ + JI^x'^ + 0{\x\^ In |x|)], (5.3) 



(5.2) 
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(5.4) 



where the relevant constants /xq, A*i, A*2 and /is are given by 

liQ = z(27r)^(aoCo - coao), A*i = z(27r)^(coai - ciOq - c?i6o), 
/U2 = (27r)^|di|^, (U3 = i(27r)^(ciai — aiCi + dibi — bidi). 

One finds the foUowing expression for the Kahler metric 

9.^ = -^ In \xf +0-^- ^J^^±l^) + 0{\x\ In \x\). (5.5) 
Then the curvature form is 

Rxx = 71 i^Tj i I . ^\oi (5-6) 
4|a;p (In |a;| + 6)^ 

where the constant C is determined to be 

I 1 2 

C = 1 _ + fti -0.738. (5.7) 

2//oAt2 2//2 

In computing Kahler covariantized derivatives with respect to ip, it is also useful to note 
that 

d^K^ = -^-^x\ii\x\'' + 0{x). (5.8) 

fJ'O fJ'O 

5.2. Distribution of flux vacua 



.0.04.. 



.0-;02- 



-0.04 



■-0.-,0-2-. 




■0.04 L 



0.04 



Fig. 3: Each point is a vacuum on the x = l — ip complex plane. The monte carlo 

simulation data is: number of random fluxes N = 5 x 10^; random flux interval 
f,hG ( — 100,100); complex structure tp space region |a:| < 0.04. There are 11249 
vacua, but 6306 of them arise at \x\ < .00001 and have been removed from the plot 
(they would all cluster at the origin). 
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Fig. 4: The plot of a numerical evaluation of 27rz ci (7/') • The monte carlo 
simulation data for each point is: number of random fluxes N = 10^; random flux 
interval f,h e (-60, 60). The data is fit by the curve - , where p = lnr. The 

conifold point r = Oisatp=— 00 for this coordinate. 
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Fig. 5: The numerical results for the number of vacua with A^aux < L for L € 
(1,972). The monte carlo simulation data is: number of random fluxes = 10^°; 
random flux interval f,h ^ (—200,200); complex structure ip space region |a;| < 
0.001. The data is fit by the curve jglooo' 
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We numerically solved for vacua around tj) = 1. To compare with the prediction from 
§2.1, we need to determine the curvature and Kahler form of the moduli space. Using the 
results of §5.1, we find (including only the leading contribution coming from the curvature 
form) 



27rz Ci('0) = const 



1 



1 



const 



(5.9) 



iM Jo (In |a;| + C)^ Inr + C 

for the portion of the fundamental region of complex structure moduli space Ai: |a;| < r 
The equations for flux vacua take the form 



DrW - 



/ • nt ^ /lap + /2^o + /3C0 
h ■ nt hiao + h2bo + ^sco 



+ 0(|a;|ln|a;|); 



(5.10) 



D^W = 



ln(a;) 



2m 
Ik 



(/i 



r/ii)(ai-^ao) + (/2-r/i2)(6i-^6o) 



/2 - r/i2 

(/3-T/l3)(ci- J^Co) + (/4-r/l4)cil 



/2 - Th2 

(5.11) 

We performed a monte carlo simulation, solving the equations of motion for randomly 
selected fluxes. In Figure 3, we plot the distribution of vacua near = 1 on the x = l — ij) 
plane. The most striking thing about Figure 3 is the clustering of vacua around the conifold 
point |]T^,0 . This would be expected from the curvature dependence of the result for the 
index (|2.1| ) (which, as shown in Figure 4, is in good agreement with our results). This is 
interesting both because it shows that the conjecture ( p.l|) contains much more information 
than a plausible boring conjecture like "vacua are uniformly distributed on the moduli 
space," and that furthermore the detailed distribution we find experimentally is in fact 
captured by the simple formula (pTTI)! 

As described in 0,0, the conifold point is particularly interesting because it pro- 
vides a natural mechanism for flux-generated hierarchies (as in the noncompact Klebanov- 
Strassler solution [0), and can provide small numbers which help with moduli stabilization 
or attaining small vacuum energy. It could therefore also be physically important that it 
dominates the distribution in Figure 3. Flux vacua near the conifold 0] (see also ||3^]) can 
explain light scalars by warping (instead of protecting their masses by supersymmetry). 
Therefore, if this feature persisted in the full problem of characterizing the string land- 
scape, it could mean that the correct stringy notion of naturalness allows warping as a 
natural explanation for light scalars (as in Randall-Sundrum scenarios |]3^), in addition to 
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low-energy super symmetry. Of course this is mere speculation at this point, but indicates 
the kinds of lessons one could hope to learn. 

As shown in Fig. 5, for the range of L most relevant to Model A of 0, we again seem 
to obtain close to scaling for the number of vacua in this region of moduli space. It 
is quite possible that at much larger L, this behavior would be modified to in accord 
with ( |2.7| ), and that a more detailed theory could explain the scaling visible here. 



6. Conclusions 

In this paper, we have seen that very simple conjectures can serve to capture the 
detailed structure of the numbers and distribution of flux vacua on Calabi-Yau moduli 
spaces. While the flux potentials are surely not the full story in string theory, they already 
encapsulate very rich behavior (including moduli stabilization, warping, supersymmetry 
breaking, etc.), and may exhibit much of the complexity of the full story. It is therefore 
heartening that simple conjectures like those of §2 can essentially classify such large en- 
sembles of vacua in simple terms. This increases our confidence that somewhere down 
the line, similar simple conjectures which capture basic aspects of the structure of the full 
space of string vacua may become available. 

The most interesting physical feature to emerge is the existence of special loci in the 
moduli space which serve as attractors, in the sense that a large fraction of the flux vacua 
cluster around them. In the one-dimensional moduli space we examined here, the conifold 
plays this role, in a way that is immediately apparent from Figure 3. The existence of 
such attractors in a more general setting could be important in helping to deflne a stringy 
notion of naturalness. 

A useful generalization of this work would involve the study of models with many 



complex structure moduli. Finding concrete models which meet the requirements of [^| 
or its generalizations should also become possible, as we gain the ability to scan over large 
.et. at oncei The e— whe.e one focu.e. on I,B flnx _ten«al. conld aI.o be 
naturally enlarged to include generic Non-Kahler type II models [|38| , or classes of vacua 
where one includes an analysis of nontrivial D-brane worldvolume superpotentials (see 
||39|| for a recent review). This could be important given that D-brane theories in flux 



^ After this paper was completed, the interesting work |37] appeared where exphcit examples 



along the lines of [13| are presented, including Kahler modulus stabilization. 
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backgrounds can give rise to very large numbers of vacua (see e.g. |5C,ri|). And finally, 
it would be valuable to study other well-motivated toy ensembles which characterize the 
sets of supergravities emerging from others limits of M-theory, like the heterotic theory. 
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